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Abstract 

In the present paper, superconvergence of second order, after an appropriate postprocessing, is 
achieved for both the two and three dimensional first order rectangular Morley elements of biharmonic 
equations. The analysis is dependent on superconvergence of second order for the consistency error 
and a corrected canonical interpolation operator, which help to establish supercloseness of second 
order for the corrected canonical interpolation. Then the final superconvergence follows a standard 
postprocessing. For first order nonconforming finite element methods of both two and three dimen¬ 
sional fourth order elliptic problems, it is the first time that full superconvergence of second order is 
obtained without an extra boundary condition imposed on exact solutions. It is also the first time 
that superconvergence is established for nonconforming finite element methods of three dimensional 
fourth order elliptic problems. Numerical results are presented to demonstrate the theoretical results. 

Keywords: Biharmonic equation; rectangular Morley element; superconvergence 


1 Introduction 

Because of significant applications in scientific and engineering computing, superconvergence analysis 
of finite element methods has become an active subject since 70’s last century. However, most of atten¬ 
tions have been paid on conforming and mixed finite element methods of second order problems, we refer 
interested readers to [21I1II1II1] for more details. Since conforming finite element methods of fourth 
order problems are very complicated, most of popularly used elements in practice are nonconforming, 
for instance, [TOl [191 [211 [25l [26] [27l [28l [29] . However, for nonconforming finite elements, due to noncon¬ 
formity of both trial and test functions, it becomes much more difficult to establish superconvergence 
properties and related asymptotic error expansions. For second order elliptic problems, there are a few 
superconvergence results on rectangular elements. In laiig, superconvergence of the gradient was ob¬ 
tained at the centers of elements for the Wilson element, which relies on the observation that the Wilson 
element space can be split into a conforming part and a nonconforming part. Due to superconvergence of 
consistency errors, superconvergence of the nonconforming rotated Qi element |20] and its variants was 
derived, see jiniiiH]. For the plate bending problem, there are only few superconvergence results for 
nonconforming finite elements. In |3], Chen first established the supercloseness of the corrected interpo¬ 
lation of the incomplete biquadratic element on uniform rectangular meshes. By using similar 

corrected interpolations as in [3], Mao et al. El first proved one and a half-order superconvergence for 
the Morley element m and the incomplete biquadratic nonconforming element on uniform rectangular 

*The first author was supported by the NSFC Projects 11271035, 91430213 and 11421101. 
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meshes. In a recent paper [5], Hu and Ma proposed a new method by using equivalence between the 
Morley element and the first order Hellan-Herrmann-Johnson element and obtained one and a half-order 
superconvergence for the Morley element on uniform mesh. That half order superconvergence can be 
improved to one order superconvergence if the third order normal derivative of exact solutions vanishes 
on the boundary of the domain under consideration. Based on the equivalence to the Stokes equations 
and a superconvergence result of Ye |80) on the Crouzeix-Raivart element, Huang et al. [7] derived 
the superconvergence for the Morley element, which was postprocessed by projecting the finite element 
solution to another finite element space on a coarser mesh. See Lin and Lin m for superconvergence of 
the Ciarlet-Raviart scheme of the biharmonic equation. Note that all of those results are only for fourth 
order problems in two dimensions. Superconvergence of nonconforming finite element methods cannot 
be found for fourth order problems in three dimensions. 

The purpose of the present paper is to analyze superconvergence of both the two-dimensional and 
three-dimensional rectangular Morley elements from |25] . Since both of them are nonconforming, one 
difficulty is to bound the consistency error. Another difficulty is from the canonical interpolation operator 
which does not admit supercloseness. To overcome the first difficulty, we use some special orthogonal 
property of the canonical interpolation operators of both the bilinear and trilinear elements when applied 
to the functions in the rectangular Morley element spaces. The other crucial observation is that the 
error between the (piecewise) gradient of functions in the discrete spaces and its mean is equal on two 
opposite edges (faces) of an element. In particular, this leads to superconvergence of second order for the 
consistency error. To deal with the second difficulty, we follow the idea from to use a correction of the 
canonical interpolation. Together with the asymptotic expansion results from [^, this yields superclose¬ 
ness of second order for such a corrected interpolation. Finally, based on the above superconvergence 
results, we follow the postprocessing idea from m to obtain a global superconvergent approximate so¬ 
lution, which converges at the second order convergence rate. It should be stressed that for first order 
nonconforming finite element methods of both two and three dimensional fourth order elliptic problems, 
it is the first time that full superconvergence of second order is obtained without an extra boundary 
condition imposed on exact solutions. It is also the first time that superconvergence is established for 
nonconforming finite element methods of three dimensional fourth order elliptic problems. 

This paper is organized as follows. In the following section, we shall present the model problem 
and the rectangular Morley element. In section 3, we analyze the superconvergence property of the 
consistency error for the two-dimensional situation. In section 4, we make a correction of the canonical 
interpolation and obtain the superconvergence result after the postprocessing. In section 5, we establish 
the superconvergence result for the three-dimensional cubic Morley element. In the last section 6, we 
present some numerical results to demonstrate our theoretical results. 


2 The model problem and the rectangular Morley element 


2.1 The model problem 

We consider the model fourth order elliptic problem: Given / G L^(r2), H C is a bounded Lipschitz 
domain. 


J A^u = /, in 12, 

( u = = 0, on dfl. 


( 2 . 1 ) 


The variational formula of problem (ED is to find u G V := iLo(12), such that 

a{u,v) := (V^M, V\)l 2 (o) = if,v)L^Q), for any vGV. (2.2) 

where denotes the Hessian matrix of the function u. 


2.2 The two-dimensional rectangular Morley element 

To consider the discretization of (12.211 by the rectangular Morley element method, let 7)i be a regular 
uniform rectangular triangulation of the domain 12. Given K G Th, let (a^i.c, a; 2 ,c) be the center of K, the 
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Figure 1: degrees of freedom 


meshsize h and affine mapping: 

4i =- 7 -> 6 =- 7 —^ , for any (xi,X 2 ) e iF. (2.3) 

n h 

On element K, the shape function space of the rectangular Morley element from m reads 


P{K) := P 2 {K) + span{a;f, xl}, (2.4) 

here and throughout this paper, Pi{K) denotes the space of polynomials of degree < I over K. The nodal 
parameters are: for any v G C^{K), 

D{v) := (^v{ai), ’ * 0 ' = 1,2,3,4, (2.5) 

where Oi are vertices of K and ej are edges with unit normal vectors Ue^ oi K, \ej\ denote measure 
of edges ej, see Figure 1. Let the reference element K he a square on (^i, ^ 2 ) plane, its vertices be 
ai(—1, —1), 02 ( 1 , —1), 03(1, 1), a4(— 1, 1), and its sides be ei = 0302, 62 = 0203, 63 = 0304, 64 = 0403. 
The nonconforming rectangular Morley element space is then defined by 

Vh := {u G L^{n) : v\k G P{K), \/K G Th, v is continuous at all internal vertices and 

vanishes at all boundary vertices, and / —— ds is continuous on internal edges 

Je 

e and vanishes on boundary edges e of Th}- 

The discrete problem of (12.21) reads: Find Uh G 14, such that 

ah{uh,Vh) ■■= (V^M/i, V^Uh)L 2 (Q) = {f,Vh)L^(a), for any Vh G Vh- (2.6) 

where the operator is the discrete counterpart of V^, which is defined element by element since the 
discrete space I 4 is nonconforming. Define a semi-norm over I 4 by 

\uh\l ■■= ah{uh,Uh), for any Uh G 14. (2.7) 


Let u and Uh be the solutions of (E3) and (ESI, respectively, by the second Strang Lemma mm) , we 
have 


\u - Uh\h S 


< C 


inf 

vheVh 


u-Vh\h+ sup 

o^WheVh 


\ah{u,Wh) - if,Wh)\ 
\Wh\h 


( 2 . 8 ) 


where the first term is the approximation error and the second one is the consistency error. Herein and 
throughout this paper, C denotes a generic positive constant which is independent of the meshsize and 
may be different at different places. 
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3 Superconvergence of the rectangular Morley element in 2D 

3.1 Superconvergence of the consistency error 

Let Ih be piecewise bilinear interpolation operator on 12, Ih ■ Vh ^ Bh, 

Ihv{P) = v{P), for any vertex P of %, (3.1) 

where 

Bh = {v€ v\k G QiiK), \/K G %}, (3.2) 

aird Qi{K) denotes the space of all polynomials which are of degree< I with respect to each variable Xi, 
over K. Let the interpolation operator he the counterpart of Ih on the reference element K. The 
bilinear interpolation opertor Ih has the following error estimate: 


\v — Ihv\Hi{K) < Ch? ^\v\h^(^k), ^ = 0 , 1 , 


(3.3) 


for any v G H'^{K). It is straightforward to see that Ih is well defined for any Wh G 14. By Green’s 
formula, 

(/, IhWh) = (A^u, IhWh) = - I "VAu- VIhWh da;ida;2. (3.4) 

Jn 


The integration by parts yields 


ah{u,Wh) 



(3.5) 


where ^ and ^ are tangential and normal derivatives along element boundaries, respectively. A com¬ 
bination of (lT4l) and (1331) yields 


ahiu, Wh) - (/, Wh) = ah{u, Wh) - (/, hwh) + (/, hwh - Wh) 

= - 'Y' / VAu-Viwh - IhWh) dxidx2 - Y] / f{wh - hwh) dxidx2 
JK JK 


Ken 


+ E 

Ken' 


f d'^u dwh 
dn? dn 


ds+ E 


Ken 
d'^u dwh 


(3.6) 


Ken 


!QK dsdn ds 


ds. 


The Cauchy-Schwarz inequality and the interpolation error estimate (IH31) lead to 



< Ch'^\\f\\L2(^n')\wh\h, 


(3.7) 


which indicates a suprconvergence rate 0{h'^). 

In the following three lemmas, we will analyze superconvergence for the three remaining terms of 

dSH). 

Lemma 3.1. Suppose that u G TLq (SI) P| i7‘*(r2) and Wh G 14. Then, 



VAu ■ V{wh - IhWh) da;idx 2 < Ch‘^\u\Hi(n)\wh\h- 


(3.8) 


Proof. On the reference element K, consider the following functional 


Bi{^,Wh)= [ d^d^, 


fK 


d^i 


(3.9) 
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Table 1: calculation of interpolation 


Wh 1 6 

6 

C16 Cl 

Cl 

c? 

cl 

ifiWh 1 ii 

6 

C1C2 1 

1 

Cl 

C2 


A simple calculation leads to the interpolations, see Table 1. 
It can be checked that 


Bl{4>,Wh) < ^1 I'/’l 

BS,Wh)=0, y^GPoiK), ywhGVh. 


The Bramble-Hilbert lemma gives 


pePo{K) 

A substitution of ^ into (13.1011 . plus a scaling argument, yield 


dAu d{wh - hwh) 


Jk dxi 

A similar argument proves 


dx\ 


da;ida ;2 < Ch‘^\u\Hi(K)\wh\h for any AT e %■ 


dAud{wh-IhWh) ^ j / ^1.21 I I I 
-dxid2;2 < Ch \u\H‘>-(K)\wh\h, 


Jk 9x2 

which completes the proof. 

Lemma 3.2. Suppose that u G i?Q (fl) p| and wt G 14. Then, 

d^u dwh 

Ken 

Proof. Given K G Th, let ei,i = I,-- - ,4 be its four edges. Define 11°. w = wds and TZ^.w = 

w — Ilg.rii, for any w G L^{K), then we have 




(3.10) 

(3.11) 

(3.12) 

□ 

(3.13) 


TTi w ds = 0. 


(3.14) 


Since ^fyds is continuous on internal edges and vanishes on boundary edges of Th, thus 


E 

KeTh' 


r d'^u dwh 
'qx dn^ dn 


= EE/ 

Ken i=i 

= EE/ 

Ken i=i 

= ztp 

Ken *=i 


d'^u dwh 
dn^ dn 


ds 


^ 0 ^ 
dn^ dn 


We first analyze the following terms 

^ J2 + J4 = 

Ken Ken 

= E 

Ken 

= E 

Ken • 


—TZ° 

dn‘^ 62 


ds • 


—TZ° 

dn^ 64 


ds 


d^u 0 dwh , 
dxj^^^ dxi . 


>€2 

pX2,c-\-h 


d^u^r, dwh , 


/ d'^u 

-h 


7T 


dwh 


dxi 


d'^u 

dxf 




dwh 


dxi 


dx2- 
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For the rectangular Morley element, we have the following crucial property 


0 dwh 


This implies 


7^ 


J 2 + J 4 — 

K&Th 


-7e° 


^0 dw. 


^0 dw. 

62 

64 

dx2 

ei 

dx2 


rx2,c+h ^ ^ ^ 

Jx 2 c-h 62 dxl ej 


KeTh 


= E 


' X2,c- 
rx2,c-\-h 


dwh 


]TZ° - 

dxi 


dx-2 




Keu 


' X2,c — h 


/d'^u 

d^u 

9wh 


62 dxj 

ej dxi 


dX2- 


The error estimate of the interpolation operators Ilg^ yields 


K&n 


K&Th 


vf — 


d^u 

] 

\dxf 

62 

~ ^ 

eJ 


LHK) 


\wh\r 


(3.15) 


< Ch^\u\Hi^n)\wh\h- 


By the same argument, we can get 


Ji + J3 < Ch'^\u\Hi(il)\wh\h- 

KdTh 


Then, a combination of (13.1511 and (13.1611 completes the proof. 
Lemma 3.3. Suppose that u S i?Q (fl) P| and Wh G Vh- Then, 

dwh 


E 

KeTh' 


dK dsdn ds 


ds < Ch \u\Hi(Q)\wh\h- 


(3.16) 

□ 

(3.17) 


The proof of this lemma can follow the similar procedure of Lemma 13.21 and herein we omit it. 
According to above lemmas, we can obtain the following error estimate. 

Theorem 3.4. Suppose that u € f G and for all Wh G Vh, then the consistency error can 

be estimated as 

ah{u,Wh) - {f,Wh) < C'/i^(||/||l2(q) + \u\Hi^n))\wh\h- 

3.2 Asymptotic expansion of the canonical interpolation 

Given AT G Tii, we define the canonical interpolation operator IVk '■ H^{K) -G P{K) by, for any 

V G 

r I I_ r 

(3.18) 


nK»(P) = v{P) and ds = d,, 


for any vertex P oi K and any edge e of K. The interpolation operator Hk has the following error 
estimates: 

\v— IIkv\h’-{k) < Ch^ Z = 0,1, 2, 3, (3.19) 

provided that v G H^{K). Then the global version n?i of the interpolation operator Hk is defined as 

IVh\K = ^K for any K G Th- (3.20) 

We need the following asymptotic expansion result from [5]. 
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Lemma 3.5. Suppose that u G H^(n), then for all Vh G Vh, we have 


h^ 


ah 


{u-UhU,Vh) < X! T 


d^u d^Vh 


3 Jx dxidx^ dxf 


dxida;2 


Ken 

sn f d^u d^Vh 
^ 3 Jii dxldx2 dxl 

^2| 


da;ida;2 


+Ch \u\Hi^Q)\vh\h- 

It is straightforward from Lemma 13.51 to derive that by the inverse inequality 


ah{u - IihU,Vh) < Ch(^ 


d^i 


dxidx2 


n(n) 


+ 


d^i 


dx\dx2 


n(Q.) 


'j\Vh\h- 


(3.21) 


(3.22) 


Based on the above analysis and Theorem l3.41 Lemma 13.51 we can get the following error estimate of 
\XlhU - Uh\h- 

Theorem 3.6. Suppose that u G H^{Vl), then we have 


\IlhU - Uh\h < Ch(^ 


d^i 


dxidx2 

Proof. It follows from Theorem 13.41 and (13.221) that 

\IlhU-Uh\l = OhiTlhU - Uh,IlhU - Uh) 


n{Q.) 


+ 




dx\dx2 


n{Q.) 


= ah{u - Uh, TlhU - Uh) + ah(JlhU - u, lihU - Uh) 

= [ah{u, XlhU - Uh) - if, lihU - Uh)] + UhilihU - u, XlhU - Uh) 
< C'^^(ll/l|L2(n) + |w|ff4(n))|II/tU — Uh\h 


+Ch( 


d^u 


dxidxl 


n(Q.) 


+ 


d^i 


dx\dx2 


n(Q.) 


|n/^u - Uh\h, 


which completes the proof. 


(3.23) 


□ 


4 Supercloseness of the correction interpolation 


In the view of Theorem 13.61 we cannot expect a higher order error estimate of \T\hU — Uh\h- To 
overcome this difficulty, we follow the idea of |3] to make a correction of the interpolation II/jU. First of 
all, we define the correction term as follows 


where Ujiv) read 


Rkv = ^ ajiv)ipj, V G H'^{K), 


(4.1) 


=-|/e. aSldx2, j = 5,7, 

“ “6 /ei dxldx2 ‘^^1’ ~ 3, 


(4.2) 


and the basis functions ipj read 


= t(a +1)^(6-1), 


= 4(6 + 1)^(6 - 1 ), 


‘^7 = t(6 + l)(a-l)^ 

2 


, ps = + 1)(6 - 1 ) 
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^2 = where we defined in (12.31) . 

Then the global version Rh is defined as 

Rhlx = Rk, for any K gTh- 

Define the correction interpolation as follows, for all v € H^{K), 

U*j^v = Ukv — Rkv, K €Th, 

Then the global version 11^ is defined as 

= foranyXer^. 

Regarding the correction term Rh^ we have the following lemma. 

Lemma 4.1. Suppose that u G then for all Vh G 14, we have 

d^Vh 
f)T^ 

KeTh ■' ^ 1 if STh ■ 

Proof. Let and ^2 be defined as in (|2.3I) . It follows from the definition of P{K) that 

d^Vh 


(4.3) 

(4.4) 

(4.5) 


a-h 


1 ■ f d^vji 1 .. f d^Vh 

{RhU,Vh) = - ^ / (05 + 07)—-3-da;idx2 + - ^ / (oe + as)-^-^-dxida;2. ( 4 . 6 ) 

jtc't; -iif 1 frc-r; -iif 2 


dx^ 


The definition of aft,(-, •) yields 


ah{RhU,Vh)= 

K&Th 


d^Vh , d^Vh ^ 1 o 


d'^RKud'^Vh f d’^Rxu d'^Vh 

dxidx2 + 2 / ——-———-— da;idx2 


(4.7) 


X dx\ dx\ 
d'^Rxu d'^Vh 


lx dxidx2 dxidx2 


dXn dXn 


da;ida;2 


(4.8) 


We are in the position to calculate three terms on the right-hand side of (14.8p . It follows the definition 
of Rx that 


f d'^Rxu d'^Vh 
lx dx\ dx\ 


da;ida;2 = 


h-^ 

~ 

h-^ 

4 

h-^ 

4 


/if L 


/if 


(6^1 -I- 2)05 -I- (6^1 — 2)07 
66 («5 + a^) + 2(05 - a^) 


d'^Vh 

dx\ 


da;idx2 


d'^Vh , ,d^Vh 


dxidx2 


dxf dx\ 

/■ d^Vh 


J d£,i{a 5 + a 7 )-^^ dxidx 2 +J 6^1 (05-I-07)11-^^ dxidx2 

h ^ /* h ^ /* 

H—^ J 2(05 - 07)-^^ da;ida:2-I- J 2(05 - ar)h-^^fi dxidx 2 . 


Since coefficients like and are constants, we can get that by parity of function and symmetry 


of domains: 


dx'f 


h-^ 

4 


/ 66(05 -I- 07) da:ida;2 = 0 , / 2(05 — 07)^16 da;ida:2 = 0 . 

Jk Jk 

ero te 

J 6^1 (as + ar)h^^^ dxidx2 j (“5 + ^ 7 )^^ dxidx2. 


Ik Jk 

Because of 05 = 07 , hence, only one nonzero term is left, which reads 


This yields 


f d'^Rxud^Vh J J 1 /■ ^ , J J 

4 = 2 4 '“=+ 


d^Vh 


(4.9) 

(4.10) 



































A similar argument proves 


d^Rxu d'^Vh 


da;ida:2 = ^ J (“6 + da;idx2. 


dxl dxl 

Note that the basis functions ^Pj{j = 5,6, 7,8) have no mixed terms, which leads to 

f O^Rku d'^Vh 


(4.11) 


^ — 0 Thus 

dxidx2 ~ 


Iji dxidx2 dxidx2 


da:idx2 = 0, 


which completes the proof. 


(4.12) 


□ 


Based on the above analysis, we can establish superclose results of the rectangular Morley element by 
the correction interpolation !!)(«. 

Theorem 4.2. Let u G Uh G Vh, be the solutions of 112.2d and i2.h\) . respectively, then we have 

|n^u — Uh\h < C'/i^(||/||l2(o) + |u|//4(Q)). (4.13) 

Remark 4.3. Comparing with the incomplete biquadratic plate element herein, the theorem does 

o3 4—1 

not require to be zero on the boundary. Besides, the correction interpolation still belongs to the 
space Vh. Because of the boundary condition ~ deduced that -i^\ = 0 and -§^1 = 0, 


ei,ej GdVL, i = 1,3, j = 2,4. Thus, = 0 and = 0, ei,ej G dLl, i = 1,3, j = 2,4. 

Proof. On the reference element K, consider the functional 


B2{u,Vh) = 


1 f d^u d^Vh 

- '-^ dCidC2 


3 Jk 

r d^i 


1 f a^u a^vh 
- '- 3 - dCid6 


'K 


a^iac. 


d6 


3 Jk acfaf2 94 
a^vh 1 

—^ d4id^2- 


a^i 


6 


'K 


a^i 


deia^2 


da 


a'^vh 

dad6 


It can be checked that 


f B 2 (u,Vh) < c\\u\\H3(K)\^h\H^K)^ 

\ B2{u,Vh)=0, yuGPsiK), yvhGVh. 

Hence, the Bramble-Hilbert lemma gives 

B2{u,Vh) < 

A scaling argument leads to 


(4.14) 


hj f a^u a^vh h"^ 

y Jk axiaxi axf ^2 + y 
,a^vh 


a^u a^vh 
ax\ax2 ax2 
a^vh 


da;ida;2 


+ ^y {a5 + a7)^g^ dxidx2 + ^ J {qq + as)^^ dxidx2 

< Ch'^\u\Hi-{K)\yh\h- 

An application of Lemma 13.51 and Lemma 14.11 yields 

ah{u - Illu,Vh) = ahiu - IlhU,Vh) + ah{RhU,Vh) < Ch^\u\H^(Q,)\vh\h- 
Then, together with Theorem 13.41 and (I4.15|l , this gives 

\X\fu-Uh\\ = ah{B.lu - Uh,Bfu - Uh) 

= ah{u - Uh, - Uh) + ah{Bfu - u, Bfu - Uh) 

= [ah{u, Bfu - Uh) - if, U*hU - Uh)] + ah(n*hU - u, Ulu - Uh) 

< C'/i^(||/||L2(n) + |u|//4(f2))|n^u - Uh\h, 

which completes the proof. 


(4.15) 


□ 
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Figure 2: macro element K 


Based on the superclose property, we can obtain the superconvergence result of the two-dimensional 
rectangular Morley element by a proper postprocessing technique. In order to attain the global super¬ 
convergence, we follow the idea of [14] to construct the postprocessing operator as follows. 

9 

We merge 9 adjacent elements into a macro element, K = |J Ki, (see Figure 2 ), such that, in the macro 

i=l 


element K, 

ywecik). (4.i6) 


We denote Zij,i,j = 1, 2, 3,4 as the vertices of the 9 adjacent elements. Then, the operator satisfies 


^3hw{Zij) = w{Zij), i,j = l,2,3,4. (4.17) 

Besides, the postprocessing operator has the following properties 

r ni,(n» = n3,u, yueH^n), 

S \'^lh'>^h\h < C\vh\h, Vuh e 14, (4.18) 

[ \u-Ulf^ulh < Ck\u\Hi(^Q), \/uGH^{Q). 


Then, we can get the following global superconvergent result. 

Theorem 4.4. Let u G Uh G 14, be the solutions of i2.2\) and respectively, then we have 

1^* - ^Ih^hlh < C'h^(||/||L2(Q) -I- |M|H4(n)). (4.19) 

Proof. It follows the properties (14.1811 and Theorem 14.21 that 


\U - Illf,Uh\h 


< |m — -I-|nfft(ni(u — u?,)|h 

< \u - Illf,^u\h + C\niu - Uh\h 

< C'/l^(||/||i;,2(Q) -I- |M|/i-4(n)), 


which completes the proof. 


(4.20) 

□ 


5 Superconvergence of the cubic Morley element 

In this section, we analyze the superconvergence property of the three-dimensional Morley element 
on cubic meshes with 17 C K.^. Let 7/t be a regular uniform cubic triangulation of the domain 17 C 
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Figure 3: degrees of freedom 


Given K £ 7/i, let 2 ^ 2 , 0 , 2 : 3 ^ 0 ) be the center of K, the meshsize h and affine mapping: 

Xl 2;i c ^ X2 X2 c ^ 2:3 X3 c r j \ r- 

4i =--» S 2 =--, 6 =--) for any ( 2 : 1 , 2:2,2:3) £ A. (5.1) 

On element K, the shape function space of the cubic Morley element reads 


P{K) := P2{K) + spanjxi, 2 ^ 2 , 2 : 3 ,Xi 2 : 22 : 3 }. 
The nodal parameters are: for any v £ C^{K), 


DM = («(<..), 


I = j = 


(5.2) 


(5.3) 


where aj are vertices of K and Fj are faces of K, see Figure 3. 

5.1 Superconvergence of the consistency error 

We also need the decomposition of the consistency error (EH). For ease of reading, we recall the 
expression as follows 

ah{u,Wh) - {f,Wh) = ah{u,Wh) - if,IhWh) + if,IhWh - Wh) 

= - ^ / VAu ■ V{wh - hwh) dxidx 2 dx 3 

Jk 


K&Th 


- ^ / fiwh - IhWh) dxidx2dx3 + 
d'^u dwfi 


d'^u dwh 


kgU-'^ 9 x 2 


ds 


+ E E 


Idxidxj dx^ 


-Tii ds. 


< Ch'^\\f\\L2(^Q)\wh\h- 


K^l~h 

A direct application of the interpolation error estimate leads to 

E / “ IhWh) dxidx2 

K(^Th'’^ 

Lemma 5.1. Suppose that u £ Hq (12) p| il‘^(n) and wt £ Vt- Then, 

/ VAu-Viwh - IhWh) dxidx 2 dx 3 <Ch^\u\ H*{n)\wh\h- 

Ken 


(5.4) 


(5.5) 
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Proof. On the reference element itT, consider the following functional 

Bs((j),Wh)= ( (5 5) 

Jk of,i 

A simple calculation leads to the interpolations, see Table 2. 


Table 2: calculation of interpolation 


Wh 1 fl 

6 

6 

C 1 C 2 

C 1 C 3 

C 2 C 3 Cl 

C 2 ' 

C 3 " c^ 

C2 

C 3 " 

IpWh 1 6 

6 

6 

Cl 6 

C 1 C 3 

C 2 C 3 1 

1 

1 Cl 

C 2 

Cs 


It follows that 


Bzijfj'Wh) c\\(j)\\j^2(^x'^\wh\fj2(^xy 

Wh) =0, G Po{K), 'iwh G Vh- 


The Bramble-Hilbert lemma gives 


B3i(l),Wh)<C inf^ \\(l} + l^\L2fK)\^h\H2fK)^^\^\m(K)\'^h\H2(^K)- 

pePo{K) 

A substitution of ^ into (15.71) . plus a scaling argument yield 

r dAudjwh - Ih Wh) 

Jii dxi dxi 

A similar argument proves 


da;idx 2 da ;3 < Ch^\u\Hi(^K)\wh\h, for any K €%■ 


(5.7) 


(5.8) 


and 


f dAud{wh - hwh) 

fx dx2 dx2 

f dAud{wh- IhWh) 


da;idx2da;3 < Ch^\u\H‘>.(K)\wh\h, 
da;idx2da;3 < Ch^\u\H‘>.{K)\wh\h, 


(5.9) 


(5.10) 

□ 


Jx dx3 dx3 
which complete the proof. 

Next, we will analyze the last two terms of (El- 

Lemma 5.2. Suppose that u G TIq (fl) P| i7‘^(r2) and wt G 14. Then, 

f d'^u dwh 

Ken i= 

Proof. Given K G Th, \et Fi,i = 1, ■ ■ ■ ,6 be its faces. Define IIp.w = wds and TZ^p.w = w — n^;..w, 

for any w G L^{K), then we have 

I TZ%wds = 0. (5.11) 

JFi 

Since /p ^^ds is continuous on internal faces Fi and vanishes on boundary faces of Fh, thus 


E x—^ f d‘^u dwfi 2i I I I 

xen i=i * * 


/ ^ ^ jrinp' fix' 
xen i=i * * 


3 6 


= E EE 


d^u dwh 


3 6 


ds 


= V W / —7?° 

2-^ 2-^2-^ Fj Q^. 

Ken i=i i=i * * 


'4 ds 
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For ease of expression, denote Ljm = Jp * = 1)2,3, j = 1, • • • ,6. 

Then we firstly analyze the following terms 


^ '^{L2 + L^)n^ 

KGTh i=l 


- E 

K€Th' 


F2 


0 dwh , , 

— dxidxs - 2_^ 


dxl dx2 


d^u Q dwh , , 

/v p, — Clx 1QX 3 


^ r^Fc+h .x^.c+h 

— h Jx3^c — h V'9^2 


F2 


Ken 
T^O 9wh 


IF 5 9x1 9x2 


F 2 


d'^u 


7^' 


9wh 


F 3 


9X2 


F 5 


da;ida;3. 


For the cubic Morley element, we have the following crucial property 


TZ 


0 9wh 


dx^ 


Fi 


= TZ 


0 9wh 

■Fi 


dxi 


0 9wh 


Fi 


7?' 

’ dx 


F2 


= TZ 


0 9wh 


dx2 


0 9wh 


Fi 


This implies 


rxi^c + h I'Xa^c+h / q2 


E E(^2 + Lb)nr = ^ 

KeThi=i Ken'^^^’‘~’^ Jx3,a-h 

_ rxi^c+h r^3,c+h 

= ^ / / 

FgT^ — Jx3^c — h 


9x1 


’ dx 


d'^u 


= rO 


F3 


dx?. 


Fe 


F 2 9x\ 


TZ 


TZ° (— 


Fi 

_ ^ 
F2 9x1 


0 9wh 
■F 2 


Fi 


9x2 

TZ 


decides 


F 2 


0 9wh 


■F 2 


9x2 


F2 


da; 1 decs- 


The error estimate of the interpolation operators 11^^ yields 

^ j2{L2 + L,)ni<ChY, I V(g; 


Ken i=i 


Ken 


v(— 


d'^u 

\ 

\ 9 xl 

F2 

9 x 1 

fJ 


< Ch?\u\H‘^{n)\wh\h- 


A similar argument proves 


L^(K) 


\Wh\} 


(5.12) 


and 


y^(Ti + L4)ni < Ch?\u\H‘t{n)\wh\f 
Ken i=i 


o 

y ] y ] (Zj 3 + Te)*** ^ Ch?\u\ H*{n)\'Wh\h- 

Ken i=i 

Then, a combination of (I5.12|l . (I5.1d|) and (15.141) completes the proof. 
Lemma 5.3. Suppose that u G Hq{^) H TZ'^i^) and Wh € Vh- Then, 


E E / <Ch^\u\H‘t(a)\wh\h- 

.^hl^JaK dxidxj dxj 


(5.13) 

(5.14) 
□ 

(5.15) 


K^Th 


The proof of this lemma can follow the similar procedure of Lemma 15.21 and herein we omit it. 
According to above lemmas, we can obtain the following error estimate. 

Theorem 5.4. Suppose that u G 7J'*(n), / G L^(n). Then it holds that 

ah{u,Wh) - {f,Wh) < C'li^(||/||L 2 (Q) + \u\H‘t(a))\wh\h, for any Wh G 14. 

We need the following asymptotic expansion result of the canonical interpolation from [^. 
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ah 


Lemma 5.5. Suppose that u G C then for all Vh G Vh, we have 

^ r d^u d^Vh 
T lr.7 

K&Th *^i=l 

It is straightforward from Lemma 15.51 to derive that by the inverse inequality 

^ d^u 


^ 1,2 p n3„, /93,, 

{u-IlhU,Vh) < ^ / a ^ 2 ~^^^i^^^^^3+Ch‘^\u\H'i{Q)\Vh\h- 


Ohiu-IlhU,Vh) < Ch X 


dxidx'p 


L2(n) 


\Vh\h- 


(5.16) 


Based on the analysis of the interpolation error and the consistency error, we can get the following error 
estimate of \U-h'a — Uh\h- 

Theorem 5.6. Suppose that u G then we have 

3 

\UhU - Uh\h < Ch X 




dh 


dxidx'^P 


L2(n) 


Proof. It follows from Theorem 15.41 and (15.161) that 
\IlhU-Uh\l = OhiTlhU - Uh,TlhU - Uh) 


= ah{u- Uh, li.hU - Uh) + a/i(n;iU - u, lihU - Uh) 

= [ah{u, XlhU - Uh) - if, TlhU - Uh)] + ahiUhU - u, UhU - Uh) 


< Ch^i\\f\\L'^(Q,) + \u\H‘i-(Q.))\^hU — Uh\h + Ch X 

i/i=l 


d'^i 


dxidx'j 


L2(a) 


\IlhU - Uh\h- 


which completes the proof. 


□ 


5.2 Supercloseness of the correction interpolation 

We can learn from Theorem 15.61 that the convergence of the error |nft,u — Uh\h is only of order 0(/i). 
Therefore, we follow the idea of [3] to make a correction of the interpolation to improve its convergence. 
The operator 11^ is modified as 

U*i^u = IIku-Rku, uGH^iK), 


_ 14 ^ 

where Rku = ^ bjiu)ipj with 
1=9 



and the basis functions 



( d\ 

+ 

d^v ') 

1 da;2da;3 

dx\dx^ 

dxidx'^ j 


f d^v 

+ 

a^v \ 

1 dsidxa 

^ dx2dxf 

dx^dx^J 

S I Fa • 

( d\ 

+ 

a^v ^ 

1 da;ida;2 

dxsdxf 

dxsdx'^ J 


' ^9 = t(a +1)^(6 -1), 
^10 = jif.2 + 1)^(6 - 1), 
^11 = + 1)^(6 - 1 ), 

^13 = j {^2 + 1)(6 - 1 )^: 
, ^14 = 1(6 + 1)(6 - 1 )^- 


(5.17) 

j = 9,12, 

j = 10,13, (5.18) 

j = 11,14, 


(5.19) 
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Cl = 


^ 6 = 


X2—X2,c 


Then the global version 11^ and Rh are defined as 


where we defined in (ED. 


RhW = Rk, for any K €%, 
n;!;|if=n^, for any if G 7 ^. 

Thus, we can ^stabilsh superclose results of the three-dimensional cubic Morley element by the correction 
interpolation 11^ u. 

Theorem 5.7. Let u G ut G 14, he the solutions of i2.2\} and i2.6]) . respectively, then it holds 

|n> - Uh\h < Ch^iWfWmn) + |u|ff4(0)). (5.20) 

Proof. On the reference element K, consider the functional 


1 


Bi{u,Vh) = '^ ^ 


^ f f f f d^u 


d^u d^Vh 

-2^7r 




K \JFi 


d^l 


d^i 


12 , 

"^JkKJf, \de^d^2 ' 563^3 

i f / f / d^u d^u 


g3 


deid6 


^3 


Vh 




dCid^d^s 




deid6 ) ^ dCid^d^. 


9CI 


It can be checked that 


J B4{u,Vh) ^ c\\u\\jj3^^~j\vh\jj2(^x)'’ 

I B4{u,Vh)=0, yueP3{K), VuheVh- 

Hence, the Bramble-Hilbert lemma gives 

B4iu,Vh) < 

A scaling argument leads to 


f d^u d^Vh 


KGTh i^j=l 


3 ./if dxidx] dxf 


1 ^ 

dxidx2da:3 + - ^ / (69-|-612)^:— 3- dxida:2dx3 

Ken 

+ 2 X! y (^10 + ^ 13 )-^^ dxida; 2 da; 3 -I-- ^ J {bn + bu)-^^ dxidx2dx3 


Ken' 

.21 


Ken 


J 

,d^Vh 

[Oio -t Oi3, 

IK ^-^2 

< C'/l"'|M|ii4(Q)|u/i|/i. 

An application of Lemma 15.51 yields 

ah{u - Illu,Vh) = ah{u - IlhU,Vh) -f ah{RhU,Vh) < C'/i^|M|ii4(Q)|z;/i|/i. 
Then, by Theorem 15.41 and (15.221) , it yields 

\lilu-uh\l = ah(n*hU - Uh,U*hU - Uh) 

= ah{u - Uh, n))u - Uh) P ah{B*j,u - u, Bln - Uh) 

= [ah{u, B*f,u - Uh) - (/, Bl^u - Uh)] + ah{Bl^u - u, Ulu - uh) 
< 0/1^(11/11^2(0) -I- |u|ii4(o))|n))u — Uh\h, 

which completes the proof. 


(5.21) 


(5.22) 


□ 
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Figure 4: macro element K 


Furthermore, based on the superclose property, we can obtain the superconvergence result of the three- 
dimensional cubic Morley element by a proper postprocessing technique. In order to attain the global 
superconvergence, we follow the idea of [14] to construct the postprocessing operator as follows. 

27 

We merge 27 adjacent elements into a macro element, K = [j Ki, (see Figure 4), such that, in the macro 

i=l 

element 77, 

nl^w&Qs{K), yw&C{K). (5.23) 

We denote Zijk,i,j,k = 1,2, 3,4 as the vertices of the 27 adjacent elements. Then, the operator Ilg^ 
satisfies 

Tllhw{Zijk) = w{Z,jk), i,j, fc = 1,2,3,4. (5.24) 

Besides, the postprocessing operator Ilg^ has the following properties 

r nijn» = ni^u, 

{ \Ul,^Vh\h < C\vh\h, (5.25) 

[ \u-Illf^u\h <Ch'^\u\H*(n), VuGi7‘‘(n). 


Then, we can get the following global superconvergent result. 

Theorem 5.8. Let u G Uh G 14, be the solutions of I12.2\} and ^2.6]) . respectively, then it holds 

\u - \ilhUh\h < Ch?{\\f\\L^Q,) + |M|ff4(n)). (5.26) 

Proof. It follows the properties (15.251) and Theorem 15.71 that 


^3h'^h\h 


< |m — -I-|n|^(n’(M 

< \u - Illi^u\h + C\iilu - Uh\h 

< Ch'^{\\f\\L^(n)+ \u\Hi{Q)), 


Uh)\h 


which completes the proof. 


(5.27) 

□ 


6 Numerical results 

In this section, we present some numerical results of the two-dimensional rectangular Morley element 
and three-dimensional cubic Morley element to demonstrate our theoretical results. Herein, we denote r 
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as the rate of convergence. For the sake of simplicity, denote 


Eti1=\u- Uh\h, ETr2=\IlhU - Uh\h, EiT3=\niu - Uh\h, 

Err4=|M-n|^n;!;u|,j, Err5=|n;l;M - Err6=|u- 

In the two-dimensional case, we choose the square domain fli = [0,1]^. We partition the domain fli 
into the uniform squares with the meshsize h = ^ for some integer N. 

• In the first example, we use the function ui{x,y) = siri^{TTx)sin'^(Try) as the exact solution of prob¬ 
lem (12.11) . 

• In the second example, we use the function U 2 {x, y) = a;^(l — x)‘^y‘^{l — y)"^ as the exact solution of 
problem (EH). 

The errors Errl, Err2, Err3, Err4 are computed on fli, the corresponding computational results of 
the two-dimensional rectangular Morley element are listed in Table 3 and Table 4, respectively. One can 
also refer to Eigure 5 for logarithmic plot of the norms above-mentioned. 

In the three-dimensional case, we choose the square domain VI 2 = [0,1]^. We partition the domain 
112 into the uniform cubic meshes with the meshsize h = ^ for some integer N. 

• In the third example, we use the function U 3 {x, y, z) = sin?'{- kx) sin?'{Try)sin?{- kz) as the exact solu¬ 
tion of problem EH- 

• In the fourth example, we use the function u/^{x, y, z) = x'^(l — x)‘^y‘^{l — y)‘^z^{l — z)^ as the exact 
solution of problem EH- 

The errors Errl, Err2, Err5, Err6 are computed on 172) the corresponding computational results of 
the three-dimensional cubic Morley element are listed in Table 5 and Table 6, respectively. One can also 
refer to Eigure 6 for logarithmic plot of the norms above-mentioned. 


Table 3: The errors of the 2-D rectangular Morley element for ui{x,y) 


N 

6 

12 

24 

48 

Errl 

3.801933642 

1.848733847 

0.916356489 

0.457125924 

r 

— 

1.040195809 

1.012556679 

1.003317321 

Err2 

1.97558386 

1.04229950 

0.526008399 

0.263522765 

r 

— 

0.922509198 

0.986612148 

0.997158239 

Err3 

1.614468104 

0.436680422 

0.111290839 

0.028096575 

r 

— 

1.886409182 

1.972243010 

1.985868661 

Err4 

3.503176938 

1.427388728 

0.359795304 

0.09010039 

r 

— 

1.295285573 

1.988130023 

1.997571101 


From the tables and figures, we can see the superconvergent behaviors of the numerical solutions. 
Besides, in our examples, the exact solution Ui{x,y),i = 1,2, or Ui{x,y,z),i = 3,4, don’t satisfy the 
boundary condition = 0,f = 1, • • • ,4, which are need for superconvergence of second order in two- 
dimensional case [iin]- However, our results still have the superconvergent property, which are coincide 
with our theoretical analysis. 
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log(errors) 


Table 4: The errors of the 2-D rectangular Morley element for U 2 {x,y) 


N 


6 


12 


24 


48 


Errl 0.014701829 0.007300451 0.003640499 0.001803668 

r — 1.009938149 1.003849379 1.013202406 

Err2 0.008054841 0.004151273 0.002093528 0.00100988 

r — 0.956302299 0.987617597 1.051752340 

Err3 0.006314902 0.001727848 0.000441448 0.00011096 

r — 1.869784037 1.968660896 1.992203814 

Err4 0.015833538 0.003874778 0.000963156 0.00024043 

r — 2.030798101 2.008272264 2.002152565 


Table 5: The errors of the 3-D rectangular Morley element for u^i^x, y, z) 


N 


6 


12 


24 


48 


Errl 4.167950479 
r — 

Err2 2.569313835 
r — 

Err5 2.221462641 
r — 

Err6 3.862477845 
r — 


1.983259265 

1.071464848 

1.299194244 

0.9837659766 

0.611765555 

1.860459095 

1.283658784 

1.589264895 


0.97524877 

1.02403111 

0.64723551 

1.00525448 

0.15663992 

1.965526946 

0.319649973 

2.005696886 


0.48523741 

1.007079492 

0.32311494 

1.002243303 

0.03948401 

1.988111511 

0.07982083 

2.001655786 



log(N) 


□ Errl 
0 Err2 
* Err3 
0 Err4 



Figure 5: The errors of the 2-D rectangular Morley element for ui{x,y) and U 2 {x,y) 
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Table 6: The errors of the 3-D rectangular Morley element for M 4 (x, y, z) 


N 

6 

12 

24 

48 

Errl 

0.001051488 

0.000509571 

0.00025198 

0.00012562 

r 

— 

1.045077305 

1.015973946 

1.004243055 

Err2 

0.000666808 

0.000330208 

0.00016451 

0.000082177 

r 

— 

1.013896 

1.00519979 

1.001368714 

Err 5 

0.00055831 

0.000154203 

0.000039494 

0.000009933 

r 

— 

1.856235564 

1.965125435 

1.991332077 

Err6 

0.00096663 

0.000232674 

0.000057562 

0.000014351 

r 

— 

2.054653764 

2.015121382 

2.003965450 


The errors of U 2 {x,y,z) 



The errors of u^(x,y,z) 



Figure 6: The errors of the 3-D rectangular Morley element for U 3 {x, y, z) and U 4 {x, y, z) 
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